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The program WTO, which is designed for computing cross sections and other relevant 
observables in the e + e~ annihilation into four fermions, is described. The various quan- 
tities are computed over both a completely inclusive experimental set-up and a realistic 
one, i.e. with cuts on the final state energies, final state angles, scattering angles and final 
state invariant masses. Initial state QED corrections are included by means of the struc- 
ture function approach while final state QCD corrections are applicable in their naive 
formulation. A gauge restoring mechanism is included according to the Fermion-Loop 
scheme. The program structure is highly modular and particular care has been devoted 
to computing efficiency and speed. 
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1 Program summary 

Title of program: WTO 

Computer: Alpha AXP-2100 (VAXstation • 4000 • 90); Installation: INFN, Sezione di 
Torino, via P.Giuria 1, 10125 Torino, Italy 

Operating system: Open VMS (VMS) 

Programming language used: FORTRAN 77 

Memory required to execute with typical data: Peak working set size: 3000-3500 

No. of bits in a word: 64(32) 

No. of lines in distributed program: 35313 

subprograms used: NAGLIB ffl 
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Keyords: e + e~ annihilation into four fermions, LEP 2, properties of the W vector boson, 
Higgs boson, initial and final state QED radiation, QCD corrections, Minimal Standard 
Model, cross sections and moments of distributions, gauge invariance, theoretical error, 
deterministic integration. 

Nature of physical problem An accurate description of the process e + e~ — > 4 fermions is 
needed in order to fully describe the physics available at LEP 2 and at higher energies. In 
particular the properties of the W boson can be correctly analyzed around and above the 
threshold only when the full gauge-invariant set of diagrams contributing to a given final 
state is included. Similarly the Z — Z production and the Higgs boson production can 
be studied within the minimal standard model. Indeed from a field theoretical point of 
view both the W and the Z bosons (and the Higgs boson too) are unstable particles and 
a full description of their production can only proceed through the complete set of matrix 
elements for the 4-fermion processes. Although at typical LEP 2 energies the difference 
between the full calculation and the double resonant approximation for a process like 
e + e~ — > / u~z7 /J z/ r r + is totally negligible, the same is certainly not true anymore when we 
consider e + e~ — > e~77 e z/ M/ u + or any of the neutral current processes. 

Method of solution The helicity amplitudes for each given process are given, according 
to the formalism of ref. 0, in terms of the 7 independent invariant which characterize 
the phase space. The phase space itself, including all realistic kinematical cuts, is also 
described in terms of invariants. Initial state QED radiation is included by means of 
the structure function approach. Upon initialization the final state QCD corrections 
are included by adopting a naive approach (NQCD). The numerical integration, with 
complete cut-availability, is performed with the help of a deterministic integration routine 
which makes use of quasi-random, deterministic number sets, the shifted Korobov sets. 
The boundaries of the phase space, with kinematical cuts, are reconstructed through a 
backwards propagation of constraints. 

Restrictions on the complexity of the problem The theoretical formulation is specifically 
worked out for massless fermions, although this is not a limitation of principle. Emission 
of photons is strictly collinear and no pt is therefore included. No interface exists with 
the standard packages for hadronization. 

Typical running time Dependent on the process and on the required accuracy. For example 
a 0.1% accuracy for a CC11 process requires something of the order of 440 CPU seconds 
for Alpha AXP-2100. On the other extreme a 0.5% accuracy for some (but not all) of the 
NC processes requires approximately 12 CPU hours on the same computer. 
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2 Long Write-up 



2.1 Introduction 

WTO is a quasi- analytical, deterministic code for computing observables related to the 
process 

e + e- -> /i7 2 /3/ 4 - (1) 

If one neglects the fermion masses there are 32 distinct processes of this kind (classified in 
ref. ||), 29 of which are at present accessible with WTO. The only exclusion is at present 
given by M f M with M = E <g> E, E = (u e , e). 

We have several codes, documented in the literature which perform a similar 

task and all of them can be classified into three broad families, i.e. semi-analytical, 
MonteCarlo (MC) integrators or classical event generator. A semi-analytical code will 
perform as much as possible of the 7(9 with initial state QED radiation included through 
structure function approach) phase space integrations leaving at most 2(3) of them for 
some kind of numerical evaluation. The advantage is represented by computational speed 
and high numerical accuracy but only very few and selected cuts can be applied, typically 
on two of the final state invariant masses. From LEP 1 experience we know that semi- 
analytical calculations are relevant also to experimentalists and not only as the ideal 
source of benchmark. A MC approach clearly represents the other extreme giving full 
cut-availability but with no intermediate analytical step taking place in the process of 
obtaining the final result out of the matrix elements. Kinematical cuts are, to a large 
extent, treated with the help of IF statements and the efficiency of the integration 
procedure is dominated by the adaptivity of the integration routine. 

WTO has been designed to lay somehow in between these two approaches. Basically 
the matrix elements are generated through the helicity formalism of ref. [[| and they 
are compact expressions completely given in terms of the invariants which describe the 
process. Also the momenta of the final states are, component by component, given in 
terms of the invariants used in the integration over the phase space, thus allowing to 
implement the kinematical cuts with an analytical control (almost completely). There 
are very few parts of the code where this is not implemented, having to do with the fact 
that in any 2^4 process there is a non-linear constraint among the invariants. Although 
this is not a limitation of principle we have nevertheless verified that with an analytical 
treatment of the non-linear constraint the balance between CPU time and efficiency in 
the integration is not in favor of the former and, for this reason, some numerical control 
of the phase space boundaries has survived in WTO. 

Having explained the quasi-analytical treatment of the calculation in WTO we spend 
now few words on the deterministic approach to the integration. The integration over the 
phase space is performed in WTO with the help of the NAG PJ routine D01GCF. This 
routine uses the Korobov-Conroy number theoretic approach with a MC error estimate 
arising from converting the number theoretic formula for the n-cube [0, l] n into a stochastic 
integration rule. This allows a 'standard error' to be estimated. There is no adaptive 
strategy at work since the routine D01GCF, being a deterministic one, will use a fixed 
grid. The evaluation of the specified observable will be repeated NRAND times to give the 
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final answer, however there is no possibility to examine the partial results but only the 
average and the resulting standard error will be printed. The error in evaluating, say a 
cross section, satisfies E < CKp~ a log a/3 p, where p =NPTS, a and C are real numbers 
depending on the convergence rate of the Fourier series, j3 is a constant depending on the 
dimensionality n of the integral and K is a constant depending on a and n. 

2.2 The phase space 

The process under consideration is specified by 

e + {p + )e-{p^) -> fMJMfsiqJJM, (2) 
and it is described in terms of 14 invariants (plus the c.m.s. energy) 



s = -(p++P-) 2 , 

xu-s = -(p + + qi- 2 ) 2 , i = 3,...,6, 

X2i-s = -{p- + qi-2) 2 , i = 3, ...,6, 

Xijs = -(qi-2 + qj~2) 2 , i^j = 3, ...,6. (3) 
Out of these we select 7 linearly independent combinations of invariants, i.e. 

m 2 _, m 2 + , M 2 , ml, m 2 , t w , t x , (4) 

such that 



x 13 


— t r — t — 

1 7 * t/ 14 W 


t l , x ls t 3 , x 16 1 t w t 3 , 


X23 


2 1 2 1 2 

= m_ + m + m - 


x 24 = 1 -m 2 -m„ -m 2 - i w +t i; 


X'25 


= m\ + M 2 + m 2 


- t 3) x 26 = -m 2 - Mg - m 2 + t w + t 3 , 


T 

•^34 


= m 2 _, x 35 = m 2 . 


, rr 36 = m 2 , 


X 45 


= Ml x 46 = l- 


2 2 71 /r2 2 2 

m_ — m + — M — m — m , 


Xc.R 


= m 2 ,. 





An auxiliary variable has been introduced, t 3 , which is fixed by the non-linear constraint. 
Additional auxiliary variables to be considered are also 



ei = m 2 _ + TOq + m 2 , e 2 = 1 + m 2 _ — m 2 + — e±, 
e 3 = m 2 + + M% + M 2 , e 4 = 1 + m\ - m 2 _ - e 3 , 

^2 = ^1) ^4 = 1 ^3" (^) 

This choice is far from being arbitrary and reflects the need of selecting those combinations 
of invariants which characterize the peak structure of the integrand. With the help of the 
previous quantities we can reconstruct completely the phase space. If Mij,Ei,9i and ipij 
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(i,j=l,4) are the final state invariant masses, the final state energies, the scattering angles 
and the final state angles respectively, then in the e + e~ c.m.s. we can write 

M 2 2 = m 2 , M 2 3 = m 2 , M 2 4 = m 2 , 

M 2 3 = M 2 , M 2 2 4 = 1 - m 2 - m\ - M 2 - m 2 - m 2 , 

M 3 2 4 = m 2 +} (7) 

and 

Ei = ^ei^/s, i = 1, ...,4, 
U 

cos 8i = 1 , i = 1, . . . , 4, 



cos -012 = 1 — 2 



2 

mi 



eie 2 

2 

m 

cos-j/'is = 1 — 2 , etc. (8) 

eie 3 

Clearly we want to take into account the QED radiation from initial states which is done 
through the formalism of the structure function |§-JK|. In this way the hard process 
becomes 

e + (x + p+)e-(x-p-) fi(qi)J 2 (<l2) f3(qs)f 4(14), (9) 
and the corresponding kernel cross section is afterwards convoluted with the standard 
structure functions. Since we want to impose kinematical cuts in the laboratory frame, 
all the quantities must be reconstructed in presence of e.m. radiation. We have 

Ei = - (x + ei - At;) y/s, A = x + - X-, 



cos 6i = 1 — 2 



X-t 



X+€i - AU 
L 7TI S 

costfj 12 = 1---^, etc. (10) 
Formally any distribution is written as 

F — J dx+dx- J dPS <d cut D(x +} s), D(x-, s) f, (11) 

where / refers to the kernel-distribution and D to the structure function. Thus / will 
depend on 

1 

X +P+ ' Qi — ~~^ X + X l,i+2S, 
1 

x_p_-gi = --x_x 2ii+2 s, 

1 . . 

Qi'Qj = --x + x_x i+2 , j+ 2S. (12) 
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In presence of radiation we write 



X — X — | — XjC — TTX . . . . . X 56 — X-^-X — TT\i | j 

X -J^2 X _ &tCm J 

x 23 = x + (ei — t 1 ), etc. (13) 



Consequently 



1 1 1 - 

X +P+ ' <?i — — x + x 13 s = — x+X-t^ = — t 1 s, etc., (14) 

2 2 2 

where s = x + X-S. f can be given in terms of the invariants {/} = m 2 _,...,t 1 if we 
normalize everything with respect to s. Indeed we obtain 

— 2 — = J dx + dx.Q cut D(x + ,s)D(x_,s)f({I},s) J dPSdj, (15) 

where 5i is the collection of delta-functions giving the invariants {/} in terms of scalar 
products. 

We now come back to the auxiliary invariant t 3 . The easiest way of obtaining the non- 
linear constraint is to start in the c.m.s of P± = x±p± with P± along the z-axis. Then q\, in 
the x — z plane, is constructed, component by component in terms of invariants. Obviously 
g 4 is fixed by energy- momentum conservation and g 2 ,g3 will be constructed recursively 
but for q 3 we allow for a fifth component. The request that this extra component be zero 
gives the non-linear constraint and therefore t 3 in terms of the remaining independent 
invariants. The non-linear constraint, being quadratic in t 3 , gives rise to two distinct 
solutions which must be properly interpreted. The phase space integral can be written as 

J dPS = J IT d% 5 + (q 2 ) 5 4 (p - £ q)j , (16) 
with P 2 = —s. After some manipulations we obtain 

J dPS = h 3/2 J U i de i U j dt j dm 2 _dm 2 U i e{e i )e{2-Y,e i ), 

x J U t d 3 qi 5(q 2 - ±e,?s) 5(X05(X 2 )5(X 12 ) J, (17) 

where 



X 1 = t, + -(ci - l)ei, 
X 2 = t 2 + i(c 2 - l)e 2 , 

X 12 = m 2 „ + l( Cl2 - l)eie 2 . (18) 
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Here c« denotes cos#j in the c.m.s. (with P + along the positive z-axis) and c\i is the 
cosine of the angle between q\ and q 2 in the same system. Moreover 



J = J d 3 q 3 5 (q% - -ep^ 5 [m 2 + -(c 13 - l)eie 3 ^ 



x 



*(l>- 



m 2 _ - m 2 + tt(c 2 3 - l)e 2 e 3 s 



To compute J we choose now another reference frame (//) given by 



(19) 



Qi = ^WI (0,0, 1,1), 

<?2 = ^2^(512,0,012,1), 

<?3 = -e 3 VI(si3COS0, Si 3 Sm0, Ci 3 ,l). 



(20) 



Let $ = S13S12 cos0 + 613612 — 1 computed for C13 = 1 — 2m 2 /e\e 3 . We get 



J 



f"K rllT 

d(j)+ d(f) 

JO J-K 



2 ei\^§ 

x 5 ^ (^2 e i - m ~ - m2 - 1 + 7^ * 

1 + fl< 



2 V5 i? 4 



(21) 



Before carrying on the full integration we can now connect the two 9 functions, #> and 
9 <) with the non- linear constraint. Let tf be the two solutions and moreover let qf(I) be 
q 3 in the c.m.s computed according to the choice t 3 = t+ or t 3 = t~ . Let finally q 3 (H) be 
q 3 as given in system 77 with 



613 = 1 " 2 



m 



eie 3 



'13 



1-c 



13- 



(22) 



To fix cos0 we first transform q 3 (I) from system I to system II and compare q\(I) with 
q 3 (II). The root for t. } is chosen to match the two vectors, one root will correspond to 
sine/) > (#>) and the other to sin0 < (#<) for fixed cos0. Moreover 



R 2 = m 2 e\ {e\e 3 — m 2 ^ s 



2 
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-e 2 (eie 3 - 2 m 2 ) c 12 + e x ^ e { - 



m 2 _ - m 2 - 1 - 2 e 2 e 3 



(23) 



where 
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eie 2 

Thus 



J dps = y J d{i}u l e(e i )e(^-J2e)je(Rl)(e > + e < ) 

x / d% 5 (q\ - ±e?a) 5 (t, + l -e x (( Cl - 1)) ff, (25) 



where {/} denotes collectively all the invariants and 



H = R; 1 I d 3 q 2 S(qi--e 2 2 sy(t 1 + t w + ^e 2 (c 2 - 1) ^ 

x5(m 2 _ + i ei e 2 (cia - 1)) • (26) 



As a final step we compute H by going back to system I. Here 



= ^ei"v/i(si,0,ci,l), 



qi 2 



q 2 = -e 1 v / 5(s2COS^,S2sin^,c 2 ,l), (27) 



and 



H = (28) 
2.3 Boundaries of the phase space 

In this section we consider the problem of reconstructing the full phase space in terms of 
our independent invariants. There are natural limits on our integration variables which can 
be derived by clustering the final state momenta qi, . . . , q± into pairs, i.e. (qi+q 2 ), (<?3 + <?4) 
or (qi + g 3 ), (q 2 + g 4 ) or (gi + g 4 ), (q 2 + g 3 ) and by imposing the phase space constraints 
relative to a 2 — > 2 process. To give an example we introduce 



We find 



P± = x±p±, Q_ = q 1 + q 2 , Q+ = q 3 + g 4 . (29) 



s = -(P + + P_) 2 , 

t = -(P + -Q-) = (m 2 _-t w ) 2 s. (30) 
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The condition to be fulfilled is that t must be in the physical region, therefore X > with 
X given by 



X = -t 2 w + (l + m 2 _ + m 2 + ) t„ 



m 2 . 



(31) 



This fixes the natural limits of the variable t w , i.e. 



T- < t w <T+, 



T 4 



1 + m 2 _ - m\ ± A 1/2 (l 



2 2 

, m_, m. 



(32) 



where we have introduced the corresponding Kallen's A-function, whose positivity in turns 
implies 



< m 2 _ < 1, < m\ < (1 -m_) 2 . 
Similarly we obtain —t x + T'_ < t. A < —t 1 + T' + with 

1 



Ti 



ei + e 3 ± A 1 / 2 ( 1, m 2 , 1 + m 2 — ei — e 2 



with the further constraint that 



A 1, m 2 , 1 + m - ei - e 2 = (e x + e 3 ) - 4 m 2 > 0, 



or -t, + T"_ < t 4 < -tj + T" + with 

1 r 



T". 



I \l/2 /i 2 2 2 2 

ei — e 4 ± A 7 1 1, ei — m_ — m , e 3 — m + — m 



with the additional condition that 



A (l, ei — m 2 _ — m 2 , e 3 — — m 2 ) > 



(33) 



(34) 



(35) 



(36) 



(37) 



The other constraints will be described directly in the presence of kinematical cuts. First 
of all cuts on the final state invariant masses are very easy to implement thanks to our 
choice of integration variables, eq.(0). Next we introduce a threshold for all energies. 
Thus if we require 



Ei > k y/s, 



(3? 



this will be the same as 



/\ ^ X-j-C^ ^ If,- 

Let us discuss explicitly the case of A positive, we get 



(39) 



X 



< 



Xe t 
A ' 



A 



(40) 
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Requiring consistency among these four conditions we obtain 



l-t w - Xe 4 + L 4 < t 3 <Xe 3 -L 3 , 
t w -Xe 2 + L 2 < t, <Xe 1 -L 1 , 
1 - X(e 3 + e 4 ) + U + U < t w < X{e x + e 2 ) ~U- L 2 , (41) 

with a derived condition on the radiated energy, 

(42) 



x+ - x _>l(l+ y £ L^j 



Next let us require that q < cos^j < Cj, where #j is the i-th scattering angle. First we 
define 

7i = ^(l-ci), r, = i(l-a). (43) 
In terms of invariants the constraints become 

e^e* < ti < Eie h (44) 



with the following definition 



J2i 1 ~~ M X + 

Eq. (|4*2])will translate into another set of constraints 



- = 1 + -^ , 

^ 1-7* x + 

i r x 

1 + — . (45) 



*3 


> 


max(cr 3 e 3 , 1 - 


£ 4 e 4 t w ) , 


*3 


< 


min (£ 3 e 3 , 1 - 


<r 4 e 4 —t w ), 


tl 


> 


max(fxiei,t w - 


- S 2 e 2 ) , 




< 


min(Ei,ei,t H/ 


- <r%e 2 ) , 


W 


> 


max (1 — E 3 e 3 


- £ 4 e 4 , ffiei + a 2 e 2 


W 


< 


min (1 - (j 3 e 3 - 


- a 4 e 4 , Siei + S 2 e 2 ) 



The optimal use of a deterministic integration routine (fixed grid) is largely based on the 
possibility of backward propagation of the constraints. Typically we deal with an integral 

J dxi 8 {pi -xx)(xx-ax) / dx n 6 [b n (x u . . . , x n -i) - x n ) 9 [x n - a n (x 1 , . . . , x n -i)] . 

;^ " (47) 

Solving the constraint b n > a n in terms of x n -\ gives c n (xi, . . . , x„_ 2 ) < i„_i < d n (xi, . . . , x„_ 2 ), 
so that we have now to solve 
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( x l> ■ ■ ■ j x n-2) < b n ^i[Xi, . . . ,X n -2) , 
C n (Xi, . . . 

On-1 (^l) ■ ■ ■ j %n-V < d n [Xi, . . . , X n -2) , 

c n (x 1 ,...,x n _ 2 ) < d n (Xi,...,X n _ 2 ), (48) 

and to iterate until the most external integration is reached. In our case all the constraints, 
including the non-linear one, are at most quadratic in the integration variables so that a 
complete solution of the procedure is available. However the rapid growth in the number 
of the constraints, expecially when several kinematical cuts are imposed, and the eventual 
occurrence of some splitting of the integrals ([a < X{ < b] U [c < Xi < d]) suggests an 
intermediate, hybrid, procedure where some of the constraints are left for a numerical 
control. Sometimes it is more convenient, in terms of computational speed, to increase 
the number of points in the fixed grid and to lose some of them instead of having a smaller 
number of points and full efficiency. This is also based on the fact that computing the 
amplitudes is usually faster than implementing the next level of backward propagation of 
the constraints. 



2.4 The Matrix Elements 

The generic process that we want to compute is characterized by six external fermionic 
legs and by a considerable number of diagrams, up to 64 diagrams for final state identical 
particles not including electrons or Higgs boson exchange. Already long time ago several 
authors JTT| realized that it is extremely unrealistic to write down the full matrix element 



M, to square it and to sum over polarizations. As a consequence of this fact various 
groups have developed their own way of evaluating the helicity amplitudes, each as a 
single complex number to be squared numerically. Our approach is essentially based on 
the request that everything has to be expressed in terms of the invariants which specify 
the process. To briefly summarize the main features of the method we first define 



m a (p) = n A u(p), u x (p) = u(p) n_ A , 
vx{p) = n_ A w(p), v\{p) = v{p)U\, 

n A = ~ (l + A7 5 ) , A = ±1, (49) 

where u, v are Dirac (massless) spinors. The key ingredient in evaluating helicity am- 
plitudes is given by the following set of formulas where all the relevant operators are 
explicitly given. 

vx(p)u a (q) = -(2 P . ? )- 1/2 n_ A [M>i + M^. 

u\{p)u a (q) = -(2p-q) '~ n A [S\,aWi + o~\ ,-o-pfl] , 

Vx(p)M0) = -(Zp-qy^Il-xlh^M + ^-vM, (50) 
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where we have introduced an auxiliary vector n M which satisfies n-n — 1, n-p — n-q — 
and it is otherwise arbitrary. In the previous relation an overall phase has been constantly 
neglected, i.e. these relations are exact modulus an arbitrary phase, for instance the phase 
of U\(p)u\(q). This fact alone will be of no consequence as long as one remembers to 
organize the calculation of all the diagrams contributing to a certain helicity amplitude 
(for a given process) in such a way that the unknown phase remains as an overall phase. 
Thus in WTO all the processes e + (p + )e~(p_) — > /i(?i) j 2(92) fzio.?) J 4(94) are organized in 
such a way that the requested operators are always relative to the pairs q 3 ), (p + , q 2 ) 
and (<7i, g 4 ). Following this strategy and given that we still want a final answer completely 
expressible in terms of invariants, we have made a very specific choice of the three auxiliary 
vectors which we need. They are 

"V = TfW ?2?4 ?1 , 

»fy = ^^nuapq^P-Qs, ( 51 ) 

where the proper normalization has been included, i.e. 

Nl — ^ { X 14 X 15 X 4b) ^ ' 
-^m = 2 ( 3 '34 a '36 a '46) > 

— 2" (Va 1 *) 1 ' 2 ' ( 52 ) 

Actually only certain combinations of the auxiliary vectors appear explicitly in the cal- 
culation. For instance $_fi</[ 3 etc. We make use of the relation 

7%^a/3 = lvlalf3l 5 + 7 5 (<S„a7/9 ~~ ^Ta + , (53) 

to obtain 



1/2 



1/2 



where again all particles are strictly massless. It is relatively easy to show that no sin- 
gularity will be introduced by the normalization factors at the boundaries of the phase 
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space. However from a numerical point of view strong cancellation may occur for par- 
ticular processes, expecially when both / and n will appear giving a factor in the 
amplitude. To avoid numerical instabilities we have made a systematic use of the freedom 
in choosing the auxiliary vectors, namely 

7 — ^ c h v ' rP-rP 

= -rre^apKP-d, (55) 

where fc a is any linear combination of P-,qi, (fe, q± , fc& of p + , p_, q 2 , (fe and k c of p + , gi, g 2 , Q4- 
Whenever a region of the phase-space is examined, which is dangerously closed to some 
portion of the boundary, then we have switched to a choice of the auxiliary vectors free 
of numerical instabilities in that region. Let us discuss this point in more details since 
it is connected with another important issue. Whenever we change auxiliary vectors for 
a particular diagram we must also remember that our amplitudes are defined modulus 
a phase, thus to implement this procedure we have to compute the relative phase. For 
instance we redefine 

N, = -^e^p^ptqi, (56) 
and assume that for a given diagram we need to switch from 



v x (q2)fv\(p+) -> vx(<h)pv\(p+), 

u_ A (p-)j4u_A(ft) -> u_x(P-)ftu_x(q3)- (57) 
First of all we introduce 

vx(Q2)MP + ) = e'+'Xt, 

vx{q2)lv x {p + ) = e i4,L X L , (58) 

and we compute 

e ^-M Xl X L = tr(^ + U 2 ), 

Xf(Xl) = tr[U)i> + UU 2 )- (59) 
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In this context it is easy to show that 



A4>i-<Pl) 



-,i{4>m—<t>M) 



-i{4>n— 4>n) 



2 ( nr T T 



(x 15 a;24 + x i5 - x 14 x 25 - 4i\e 4 ) , 



1/2 (^34 X 56 X 35 X 46 X 36 X 45 ^A^s) , 



2 (t np np rp \^ I 

V J '34 a '35 a '46 a '56 / 

£15 = e(5l,5i2,53,54)- 

A typical example of a diagram to evaluated will be 

D(X,p,a) = U A (p + )rM_ A (p_)M p (gi)r / ^ p (g 2 )M C (g 3 )r"t;_ 5 (g4), 
where T, . . . T" are strings of 7— functions. This will be computed as 



(60) 



(61) 



D(X,p,a) = tr {V [v- p (q 2 )vx(p+)]T [u_ A (p_)^(g 3 )] T" b- ? (54)«p(5i)]} 



= E 

A,p,S=-l 



i=l,5 



(62) 



The presence of saturated e-tensors seems to represent an obstacle towards giving the full 
amplitude in terms of invariants. Our strategy therefore is to express all \e\ in terms of 
invariants. For instance if we define S\ = e(p + ,p~, q±, q 2 ) then we obtain 



I £> I 2 O / ~j ,y» /->-» /->-» I rv* rv* rv* I /v» ry* rv* \ 

|C| — ^ \ J '13' l 14 J '23 J '24 T i 13 J '24 J '34 T J, 14 J '23 -I '34/ 1 



2 2 2 2 2 

13 24 14 23 34 * 



(63) 



and similar expressions for the other 5 independent saturated e-tensors. 



£1 = £(p+,P-,5i,52), 

e 2 = e(p+,p_,gi,g 3 ), 

£3 = £(p+,P-,52,53), 

£4 = £(p+,5i,52,5s), 

£ 5 = e(p_, 51,52,53)- 



(64) 



Moreover the sign of £162, • • • , £i£5 can all be determined in terms of invariants. Let rj be 
the unknown sign of £1, then we end up with some expression of the form 



a kp£ b \,P,a si s n ( £ i^) \ £ i 

i=l,5 



(65) 



Therefore, as long as we only need |<i| 2 , there will be no problem at all with the unknown 
r] since rf = 1. Of course we can have additional imaginary parts in the amplitude, due to 
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the width of the vector bosons. Suppose that in a diagram some internal line can be either 
a photon or a Z. Schematically we denote with a v + ib v e, V = 7, Z the corresponding 
amplitude and with a + if3 the Z propagator. The overall expression for the amplitude is 
therefore 

d = (a 7 + aa z ) — f3b z e + i (5 7 + I3a z ) e + i f3a z1 (66) 

In computing |<i| 2 all terms which are linear in e do not contribute after integration, 
therefore 



\d\' 2 = (a 7 + aa z ) 2 + (/3b z e) z + (6 7 + (3a z ) 2 e 2 + {(3a z ) 



(67) 



We give an example of how the helicity amplitudes are worked out explicitly. Let us start 
from some diagram whose expression is 



x u P (^h U l+ (h +$-- Aa) 7 M (a + &7 5 ) v- p (q 4 ) 

x ^(gih^+^te), 

where 7 ± = 1 ± 7 5 . It is easy to show that 



(68) 



iq (A, A, A) = 16 g 2 L 5 x ,i (x 15 x 2 aX3aXa5Xa6) 



uj (A, —A, —A) 



where 



-1/2 



x tr 



n_A7ViU2^ 0+ + *-- Aa) ri%&2&ii u i%& 3 ^ 



= 16 9 5x ,1 (x 15 X2AX3AXA5 x A6 



,-1/2 



x tr 



+ 



oj(X,X,-X) = oj(X,-X,X) = 0, 



^>9 L 9 R h,-l (^14^25^34^46) 



-1/2 



x tr 



n-AT^-jkV (i4 + + 14_ - Aa) i^AAA^AA 



= l6g L 9j\-l(x 14 X 2 5X 3 AXA6) V2 



x tr 



iUt^t" + i>- - Aa) i^AaAA^A^ 



9l = \ ( a + b) , g R = ^(a-b), 



(69) 



(70) 



after which the Kahane identities can be applied and the trace can be taken. For instance 
we find 



lo (A, -A, -A) = -32 g L g R 5 x ,-i (xu^sXuXw) 



'1/2 



x tr 



n_ x ri>_A 3 AA2<hr 0+ + i>- - Aa) A 2 i> + 



= 64g L g R 5x~i(xiAX 25 x 3 AXA6) 



-1/2 



x tr 



^-xAaAiAM- {i>+ +i>-~ Aa) 



(71) 
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There is another place where we must worry about relative phases. Consider the NC64 
processes, typically e + e~ — > / ///. The 64 diagrams arise from 4 permutations of 16 
different topologies, which we denote by 

92,93,54} = {91,92,93,94}, 

^{91,92,93,94} = {93,92,91,94}, 

^3{9i,92,93,94} = {91,94,93,92}, 

^i{9i,92,93,94} = {93,94,91,92}- (72) 

In the massless case there are 12 distinct sets of helicity states which we denote by labels 
h = a — m. For instance the effect of V2 can be inferred from table 1. 



p+ 


V- 


9i 


92 


93 


94 


h 


V 2 ^ 


P+ 


V- 


9i 


92 


93 


94 


h 


+ 




+ 




+ 




a 




+ 




+ 




+ 




a 




+ 




+ 




+ 


b 






+ 




+ 




+ 


b 


+ 






+ 




+ 


c 




+ 






+ 




+ 


c 




+ 


+ 




+ 




d 






+ 


+ 




+ 




d 


+ 




+ 






+ 


e 




+ 








+ 


+ 


i 




+ 




+ 


+ 




f 






+ 


+ 


+ 






j 


+ 






+ 


+ 




g 




+ 




+ 


+ 






1 




+ 


+ 






+ 


h 






+ 






+ 


+ 


m 



Table 1: The effect of the P 2 permutation on the non-zero helicity sets for the 16 funda- 
mental topologies contributing to a typical NC64 process. 



Our strategy has been to compute explicitly the 16 diagrams corresponding to the stan- 
dard ordering {91,92,93,94} and to perform the relative permutations of invariants and 
helicity sets for the remaining diagrams. This however requires the knowledge of all the 
relative phases which we have computed explicitly. Since the relative phases depend on 
the particular assignment of momenta and helicities but not on the particular topology 
under considerations it is enough to compute them for one diagram per permutation. 

As we have already stressed all fermions in WTO are kept rigorously massless. This 
fact however should not be confused with a limitation of the helicity formalism. Indeed 
we briefly sketch how fermion masses can be incorporated into our approach. We con- 
sider spinors u(p, n, A), A = ±1 defined to be the eigenstates of the operator P+(p, n, A) 
corresponding to eigenvalue 1. 

P+(p,n,X) = A + (p)^(l + iA 7 V)A+(p), 

A+(p) = ^-(-^ + m) 7 4 , (73) 

where n is the polarization vector and p ■ n = 0, n 2 = 1. If we consider the typical CC20 
process, e + e~ — > e~I7 e z/ M/ u + and we are interested in the finite m e effects in the limit of 
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zero scattering angle then things can be organized in such a way that we only need to 
compute 



u (p, n, a) u (q, A) , etc., (74) 
where p 2 = — m 2 while q 2 = 0. In this case we easily obtain 



u(p,n,a)u(q,X) = \{mp q ) 1/2 



8 



]T N a (p, q) H t + m) (1 + ia i) i (l + A 7 5 ) V a (X, a), 



x 

a=± 

mN~ 2 = amq-n-p-k, V a {\, a) = - (1 + aXa) . (75) 

2 

Moreover if both fermionic lines are massive we can repeat the whole procedure by using 

u(p,n,-X)u(q,n,X) = (-ifi + m) (l - iX^ii) + m) (l + ?A7 5 )d) , 
u(p,n, X)u(q,n, A) = --^ (-■ ip 1 + m) (^1 + iA7 5 )d) ^ (-^ + m) (l + iX^ 5 ^ , (76) 

with 

N 2 = (p + qf. (77) 

For NC32 and non-leptonic NC64 or Mix43 the exchange of gluons must be taken into 
account. This raises the basic question of what to use for a s or, alternatively, at which 
scale to compute it. Here the situation is far from satisfactory, expecially as compared 
to the well established predictions for 2-fermion production at the Z resonance. In the 
present version of WTO we have adopted the ad hoc prescription of using a s (M w ) which 
becomes therefore an input parameter. However a more satisfactory choice would imply 
to consider a running a s . Ideally for each sub-process g*(s) — > qq we should include a 
factor a s (s) and subsequently integrate over the invariant mass s, however this can only be 
done when a reasonable cut is applied, s > s to avoid entering non-perturbative regions 
where the correction factor would go out of control. Clearly more refined analyses are 
needed to reach a satisfactory level of description. We end this section with a very short 
description of what it is meant by NQCD. Consider the CC11 process e + e~ — > fi'u^ud 
which effectively is a CC10 process. One would like to include final state QCD corrections, 
even when kinematical cuts are imposed. By naive QCD we mean a simple recipe where 
the total H^-width is corrected by a factor 

r„ r„, (i + 1 «m , (78) 



and where the cross section gets multiplied by 



/ a s (M w ) \ 

&ccw,c — > o- C cw,c I 1 H I • (79) 
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This naive approach, consequence of our ignorance about the complete result, would 
be correct only for acco3,ex, the double- resonant approximation with fully extrapolated 
setup. For accio,c it is instead only a rough approximation because of two reasons. First 
of all in CC10 we have not only a virtual QCD correction to the Wud vertex but also a 
box diagram. Moreover QED and QCD radiation are quite different if cuts are imposed, 
expecially in presence of severe cuts. Thus any inclusion of final state QCD corrections 
is, at present, only a very crude approximation which moreover can become quite bad 
whenever stringent kinematical cuts are applied to the process. 



2.5 Mappings of the phase space 

For a total cross section we need a 7-dimensional (9-dimensional when initial state ra- 
diation is included) integration. In order to achieve the requested precision we must 
examine the structure of the integrand which will show a complicated peaking structure 
due to propagators (both in the time-like and in the space-like region) , to Jacobians and 
to square-integrability but not integrability of the structure functions. As a consequence 
we must perform a careful analysis and find out, as much as possible, those isomorphisms 
of the phase space which are needed to cure the peaks of the differential distribution. 
First of all by transforming the integration over the phase space into an integration over 
the kinematically independent invariants we have introduced a Jacobian which contains 
two inverse square roots of polynomials. This will always be cured as follows. The part 
of the Jacobian denoted by R 2 is a quadratic form in t 1 so that our cross section can be 
schematically written as 



a 



dm 



dt 



ft 

J dt 1 f(...,t 1 )(-at 2 + 2bt 1 + 



-1/2 



(80) 



where a,b,c are functions of the other invariants and moreover a > over all the phase 
space. We introduce the following mapping 



1 b — at 

t = ■= arccos 

'a D 



D 



b + ac. 



which gives 



dm 2 _ 
1 



dt 



J dt\ f i^. . . ,b/a + D/asi^y/ati 



arcsm 



a 



arcsm 



a 

b D 

~T — • 

a a 



l + -(t+- U) 



1 + ^ 



(82) 



There is another square root in the Jacobian which is slightly more difficult to treat since 
it is a polynomial of order > 2 in all invariants. Actually it can be seen as a quartic form 
in the variable z = m 2 defined in eq.(H), i.e. 
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R\ = - [z 2 + (E - 1) z + A+] [z 2 + (E — l)z + A_] , 
where we have introduced 



E 



Given 



r l,2 
Si, 2 



m 2 _ + m 2 + + Ml + ml, 
M m ± m_m + . 



2 I, 1 " S =F V A + ) » 
(l-E) 2 -4Ai, 



we have the following three possibilities 

A_ < 0, A_ > > A+, A + > 0, 

which correspond to 



K < 0, 

T\ < z < Si, or s 2 < -2 < r 2 . 
Only the second one will be discussed here in some details. Again schematically 



n > ri, r u <r 2 , r lj2 = r =F - yA. 



Next we introduce a new variable 



cos = 2 



r — z 
VAT' 



and distinguish among three alternatives. If r u < r then 

I = J dfjLf\r- -yA_cos0, . 
If instead r\ > r then 

/•w / 1 / 

1= dfif r + -JA_cos0, . 
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fji = -F(<f>,K). 

7 



fi = -F((/>,K). 

7 



Finally if T\ < r < r u we first introduce 



and subsequently we obtain 

J = jT d///(r-- V / Alcos0,...) + jT d^f(r + ^^A. 



COS( 



In the last equations F denotes the elliptic function 



and 



F (0, K) = J # (l - K 2 sin 2 if) 



-1/2 



(92) 



(93) 



(94) 



R 2 — 


-(z - r 1 )(z - r 2 ) 


'(* 


A = 


s 2 + (r 2 - r)(ri - 


r) 








K x = 


A+(A 2 + l) 1/2 , 




K 2 = 


1 








5 






7 


/ ^ \l/2" 

(2sk 1v ^t) 



(95) 



The case of four real roots will not be considered explicitly here but it can be treated 
along the same lines, i.e. by introducing elliptic functions. 

Whenever initial state radiation is included we have to consider the following two 
integrations 



a = J dx + dx- D(x + , s), D(x_, s) J dPSQ cut f, 



where the structure functions are 



D(x,s) = '-(l-x) 



P „ _^/ 2 -i I ex P 



¥ (! - i E ) 



r (i + \p) 



D r (x, s) 



(96) 



(97) 



and the /3-factor is 



n \ mi 



(98) 



In order to have a proper treatment of the x + /x^ integration we first change variables to 

v 



U 



(99) 
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The behavior near x — 1 is cured by introducing two new variables, < [/, V < 1 which 
are related to u, v through 

u = 1 - (1 - «,) (1 - £/) 2//J , 

v = u - (u-u,)(l -\/) 2//3 , (100) 

where we assume that the final state invariant masses are bounded by below, i.e. m 2 : > 
Mj 2 , . . . ,m 2 + > Mf and 



^ = max(u a , u b , u c ) , 

u a = (Mi + Mg) 2 , 

u b = (M 2 + M 5 ) 2 , 

Mc = (M3 + M4) 2 . (101) 

Finally most of the diagrams belonging to a given process receive resonating contributions 
by vector boson propagators. The resonating mapping is not always mandatory, expecially 
for those processes where the final answer is not dominated by W or Z peaks, like for 
instance in NC processes where we can have large effects from photon exchange or from 
t-channel diagrams. However for CC processes we have a double-resonating H^-exchange 
which will show up as 



a = j . . . j dm 2 _ [ dm 2 ^ 



X 



.../(m 2 ,m+...) 



(vm 2 _ - ^) + {ym_ lw f (vm 2 + - /i^) + {vm +lw ) 

(102) 



where n 2 w = M^/s and ^ w = T w /M w . The lower and upper limits of integrations for 
m± may differ from the natural limits because of kinematical cuts. Let us deal with one 
specific integration and change variable, e.g. z = vm 2 _. The typical integral will be 

/(*,.••) 



1 r z + 
I=- dz 

J z 



We first transform 



/4 



1 + ±W- 
1 ^ 1 M,„ 



V J*- (z-fl 2 ) 2 + (Z lw f 



1 + jf- tan [(z b -z a )Z + z a ] 



(103) 



(104) 



with 



z a> b = arctan ■ 



w 1 \ M, 



(105) 
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and where < Z < 1. In this way we get 

I = ^-Za j l dZf{z{Z) ^ ^ (1Q6) 

Concerning the resonating mapping we stress that indeed there are 4-fermion processes 
which receive the largest contribution from double-resonating diagrams. A typical ex- 
ample is given by e + e~ — > yTV^udd at LEP 2 energies where the percentage difference 
between CC03 and CC11 is totally negligible. In this and similar situations the isomor- 
phism of the phase space that we have described is mandatory. The other cases where 
the mapping is strongly advised are those where very stringent cuts are applied around 
the double resonance. For instance 



e+e" -> bbff, M z - A < M(ff) < M z + A, M(bb) > M . (107) 

On the contrary there are examples where double-resonating diagrams only give a small 
fraction of the total answer, either because there are 7-exchanges not suppressed by 
kinematical cuts or because of large t-channel contributions. In all these cases it turns 
out to be more convenient to deal with a flat phase-space, i.e. no mapping at all. 



2.6 Gauge invar iance 

Any calculation for e + e~ — *> 4-fermions is only nominally a tree level approximation be- 
cause of the presence of charged and neutral, unstable vector bosons and of their inter- 
action with photons. Unstable particles require a special care and their propagators, in 
some channels, must necessarily include an imaginary part or in other words the cor- 
responding 5 1 - matrix elements will show poles shifted into the complex plane. In any 
field-theoretical approach these imaginary parts are obtained by performing the proper 
Dyson resummation of the relative two-point functions, which at certain thresholds will 
develop the requested imaginary component. It is well known that the bosonic parts of 
vector boson self-energies are not separately gauge independent, or stated differently they 
will show a gauge-parameter dependence. This problem will arise in the context of a 
full one loop calculation, while here we can take into account only the fermionic parts 
which are obviously free of ambiguities. Thus the correct recipe seems representable by a 
Dyson resummation of fermionic self-energies where only the imaginary parts are actually 
included. As a result the vector boson propagators will be inserted into the corresponding 
tree level amplitudes with a p 2 -dependent width. Its has already been noticed by several 
authors that even this simple idea gives rise to a series of inconsistencies, which some- 
times may give results completely inconsistent even from a numerical point. The fact is 
that the introduction of a width into the propagators will inevitably result, in some cases, 
into a breakdown of the relevant Ward identities of the theory with a consequent viola- 
tion of some well understood cancellation mechanism. Sometimes the effect of spoiling a 
cancellation among diagrams can result into a numerical catastrophe. 

In the present version WTO is including reparation of gauge invariance according to 
the fermion-loop scheme [0 for the CC20 family, while an extension to other processes 
in the TeV region is currently in preparation. 
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Here we briefly sketch the procedure adopted for CC20. In this case there is a violation 
of gauge invariance induced by the so called fusion diagram which contributes to e + e~ — > 
e~V e ud through a "f(Z) and W bremsstrahlung followed by 'y(Z) + W — > W* — > wrf. 
In WTO the width of a vector boson is always zero for a t-channel exchange thus for 
CC20 gauge invariance will be violated because of the W width in the s-channel. The 
fermion-loop scheme includes in the calculation the imaginary part of the two triangle 
diagrams (with opposite charge flow) with a fermionic loop which represents the first 
order correction to the VWW vertex. In the fermionic loop we include leptons, neutrinos 
and u, d, c, s quarks, i.e. everything which is allowed in the W decay. In principle and for 
m 2 + > (mt+rrib) 2 there are two additional imaginary parts for the vertex which correspond 
to cutting a t-line and a 6-line. These extra, m^-dependent, contributions have not been 
included so far in our calculation. To be more specific we take into account the corrections 
to 



V v (q) + W ll { P ) - W a (k), V = 1 ,Z 

V = P--5i, P = P+-Q2, k = q 3 + q 4 . 



The correction factor becomes 



16 IT 



(108) 



(109) 



where we have explicitly factorized the lowest order coupling. The corrected vertex V does 
not factorize into the corresponding lowest order amplitude but instead three additional 
form factors must be included. They give 



fj.ua 

Q 

c 4 



=1,4 



y 

A 

y 

A 



di + y+ j -^0 + \ Ci + y j L 



i = 1,...3 



' b 4 eA ( d 4 fA 



A A 2 



A A 2 



(110) 



with 



r2 



^3 



Moreover 



Vuan = "^^^f ( U1 ) 



24 



A = (p.qY 
h = 



2 2 

P 9 , 



log 



p - q 



-Vx 



(112) 



and 



4 



ai 


= 


6x - -z, 


a 2 


= 




a 3 




3 


/7 A 




8 


6i 




—x(yz — Axy + xz), 






3 / 

— xl — ?vz 4- 4t?/ — tzi 


h 


= 


■y 

-x(5yz — Ay 2 — 8xy + xz), 






1/176 ,„ 9 4 , 






2 J/( 3 ^y + 16xz I6x + yz), 


Cl 




-4x(y - z + 2x), 


c 2 




2x(-y + 2; - 4x), 


c 3 




2a;(3y-2 + 2a;), 


c 4 




-16xy, 






— x 2 (— 3yz + 2y 2 + 6xy — xz), 






3 

-x 2 {3yz — 2y 2 — 6xy + xz), 
2 






-x(—2y z — 9xyz + \Axy + 10:r y — x z), 






-4xy(-3yz + 2y 2 + 18xy - 3xz + 2x 2 ), 


e 4 




—5xy 2 (yz + 3xz — 6xy — 2x 2 ), 


/4 




-10x 2 y 2 (-2yz + y 2 - 2xz + Qxy + x 2 ), 



with 

x = q 2 , y = k 2 , z = k 2 + q 2 -p 2 . (114) 

We have verified that the inclusion of the full correction in WTO accounts for an increase 
of the CPU time corresponding to a factor of 1.3, which we consider still very reasonable. 
It should be stressed here that we include the imaginary part for both ^WW and ZWW 
vertices in order to preserve the full SU {2)®U (1) gauge invariance. This can be most easily 
seen by considering the sub-process e~W + — > e~ud and by writing the corresponding 
Ward identity p^A*^ + i^ 2 M w A ( p = in the gauge (0 being the Higgs ghost). The 
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relevance of the fermion-loop scheme for CC processes with electrons in the final state 
becomes more and more pronounced for smaller and smaller electron scattering angles. 
We have explicitly verified that a selection cut of 10° on the e~ scattering angle in the 
process e + e~ — > e _ z/ e z/ M /i + is enough to avoid any problem. Partial results are shown in 
table 2, where 6> m Z(e~,e~) and we have defined 5 FL = 1 — <j/cr FL . 



d m (deg) 


o (nb) 


a FL (nb) 




10° 


0.20140(4) 


0.20143(4) 


+0.01 


5° 


0.20653(3) 


0.20643(3) 


-0.05 


1° 


0.21554(21) 


0.21014(3) 


-2.57 



Table 2: Inclusion of the gauge restoring terms at 190 GeV in e + e — ► e v e v^n + , for 
dm < e - < 7i — 9 m . In both cases the W width is running. Moreover E e ^ > 1 GeV and 
0(e-,fi + ) > 5°. 



Needless to say the fully extrapolated cross section for e + e~ — ► e^u e u^ + {ud) requires 
both gauge reparation and non zero electron mass. Thus strictly speaking one cannot set 
9 m to zero in the present version of the program. For a better understanding of the effect 
of restoring gauge invariance in CC20 we have shown in Fig. 1 the angular distribution 
da jd cos 9i where Q% is defined by 

ti=z+ei ■ 1 + cos6 i ( 115 ) 

X- + x + — [X- — X + ) COS Ui 

i.e. 9i is the angle between the incoming e~ and the outgoing l~ . In Fig. 1 we have 
shown the angular distribution for e + e~ — ► e~V e P ll jjf r with and without the fermion 
loop corrections as compared with the similar distribution for a CC11 process e + e~ — > 



2.7 Final state QED radiation 

In this rather short section we would like to give an idea of the problematic connected with 
QED radiation. Because of gauge invariance there is no meaningful splitting of the QED 
corrections between initial state and final state, splitting which was instead applicable for 
LEP 1 physics. There are various attempts to circumvent the problem, most noticeably 
the current-splitting of GENTLE |13[1 , but as a matter of principle even this technique 
cannot guarantee an unambiguous answer. 

From this point of view only the full 0(a) calculation makes sense and correctly 
reproduces the 0(ax const) terms of the QED corrections. As in any other calculation 
of this type we would end up with a correction Sqed factor of the following structure 

0~QED = $soft + $virt + $hard, (H6) 

after which the proper exponentiation can be performed. What we know of Sqed are the 
universal, leading logarithmic terms but a calculation of the virtual corrections, including 
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up to pentagon and hexagon diagrams, is missing. In such a situation it is extremely risky 
to partially account for final state QED radiation, even though any realistic estimate of 
cross sections and of energy losses should include at least some reasonable guess. 

To this end we have preliminary investigated the inclusion of QED final state radiation 
to CC11 processes. We have made use of our knowledge of the universal soft terms and 



moreover of the general and simple result of ref. |I4| for the emission of hard collinear 
radiation by charged particles. The latter allows to include hard photons emitted within 
a cone of fixed half-opening angle 5 C , in the limit 5 C << 1. This result is based on a 
rigorous gauge invariant procedure at 5 C << 1 so that our correction factor for each 
charged fermionic line will be 



exp 



"a 

— 

7T 



soft 



1 



01 A 
— O co ll 

7T 



(117) 



In this way we can only make a very rough estimate of the effect of final state radiation 
since all the hard constant from virtual corrections are still missing. Moreover the proce- 
dure of exponentiation is also far from unique since we could decide to exponentiate not 
only the leading logarithms but also numerically relevant terms from 5 co u . Let Ej , be 
the energy and the mass of the emitting fermion, moreover we denote by E = e ^fs the 
energy threshold, i.e E > E . In our formulation 



Ei = -rjiy/s, rji = x+ti + (x+ - x_) t { . 



'118) 



According to ref. [14 and to the subsequent generalization |L5| we define 



Pi 

1 - Si 



0~cVi 

2rrii 



(119) 



and compute 5 co u as Q\ C with C given by eqs.(5-6) of ref. [|I5 



'u^ud, 



Our preliminary results show very little effect on the cross section for e + e~ — > n~ 
for instance we obtain cr = 0.59192(4) nb and o~ FSR = 0.59165(4) nb at y/s = 190 GeV and 
for 5 C = 5°. The effect of QED final state radiation must certainly be included for any 
reliable determination of the physical observables at LEP 2 but we have decided for not 
including it in the present version of WTO since a more detailed theoretical investigation 
is needed. 



2.8 Program structure 

We start this section with a brief description of the general features of the program. 
With WTO it is possible to access 29 out of 32 4-fermion (massless) processes. The most 
relevant INPUT parameters are given in the BLOCK DATA INIT. They are: 

1. WM. The W mass (GeV). 

2. ZM. The Z mass (GeV). 
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3. ZG. The Z total width (GeV). 



Other quantities like the W width are not external input parameters but rather they 
are computed internally as given by the standard model. There is no possibility of chang- 
ing the weak mixing angle independently from the others input parameters, thus s 2 w is 
always a derived quantity In addition to the 29 processes which include gluon exchange 
in the non-leptonic case there are 4 processes (e + e~ — > ffbb) which can be computed and 
where the Higgs boson exchange is included. All processes receive and internal number 
according to the scheme given in appendix. Moreover each process belong to a particular 
class, but it is important to know only the highest class in a particular family, thus CC11, 
CC20, NC21, NC24, NC32, Mix43, NC48, NC64. WTO will always perform a call to the 
NAG routine X02AJF which returns the machine precision (ZRM). Some of the internal 
controls are based on this quantity so one should be aware of the fact that, for instance, on 
a VAXstation • 4000 • 90 ZRM = O.UE - 16 while on a Alpha AXP-2100 ZRM = 0.1 IE - 15. 
It is perhaps appropriate to give here a brief summary of the computational speed versus 
requested precision for some of the processes. In table 3 we have indicated both entries 
for a process representative of each family. 



Family 


Process 


Rel. error 


CPU time 


ecu 


All 


< 0.01% 


3 : 


38 : 


53 


CC20 


All 


< 0.02% 


2 : 


42 : 


57 


NC24 




0.25% 


11 


: 47 


: 47 


NC32 


UUCC 


0.08% 


12 


: 30 


: 52 


NC19 




0.48% 


12 


: 06 


: 14 


Mix43 




0.08% 


7 : 


59 : 


26 


NC48 


e + e~uu 


0.17% 


2 : 


10 : 


21 


NC64 




0.03% 


5 : 


24 : 


42 



Table 3: CPU time needed on a Alpha AXP-2100 computer for a given relative error (in 
percent) and for the some representative processes in each family. 

As it can be seen from table 3 there are considerable differences in CPU time between 
different processes. The reason is not only connected with the number of diagrams con- 
tributing but also with the number of helicity sets required by each process. Always 
working in the massless fermion limit we have only 2(3) sets of non zero helicity states 
for CC11(CC20) but up to 8(12) for NC24(NC48). 

2.9 Input 

The meaning of the input parameters is the following. 
0PRT( CHARACTER* 1) 

There is the possibility of printing some additional information about the calls to the 
various subroutines. This however requires a detailed knowledge of the internal structure 
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of the calculation. Thus, by default, 0PRT='N' and it can be set to 'Y' only by changing 
the corresponding line in BLOCK DATA INIT. 

NPTS ( INTEGER*4) 

The actual number of points for the integration. Built-in choices are NPTS=1 — 10. If more 
integration points are needed then the array VK(NDIM) must be re-initialized through a 
call to the NAG routine D01GZF. The optimal coefficients for p-point integration over 
the n-cube [0, l] n require that p is a prime number or p is a product of 2 primes, P2 and 
Pi chosen to be the nearest prime integer to p\. The built-in choices are shown in table 4. 



NPTS 


p = pixp 2 


1 


2129 


2 


5003 


3 


10007 


4 


20011 


5 


40009 


6 


80021 


7 


99991 


8 


10193 x 101 


9 


22807 x 151 


10 


32771 x 181 



Table 4: Built-in choices for p, needed for p -point integration over the n-cube [0, l] n . 



NRAND ( I NTEGER* 4 ) 

The number of times where the integral is computed in order to give an estimate of 
the standard error (usually NRAND < 6). A 'fast' estimate of the result with standard 
error around 1% is already achieved with NPTS=4, NRAND=4. A more accurate but still 
intermediate answer will require (8, 6) and a very precise, but also very time consuming, 
estimate will proceed with (10, 6). There is no intrinsic limit to NRAND while NPTS < 10 <, 
unless an independent call to NAG routine D01GZF is performed. 

IPR(INTEGER*4) 

The catalog number of the process [1 — 33]. The internal structure of WTO is organized 
in such a way that the user does not have to bother about quantum numbers of the final 
state. 

IPR0(INTEGER*4) 

If IPR < 3 there is the additional possibility of taking into account only the CC03 part 
of the CC11 family. 

RS(REAL*8) 

The c.m.s energy in GeV. 
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OPEAK (CHARACTER* 1) 

In general we are dealing with double resonant, single resonant and non-resonant diagrams. 
0PEAK='Y[N]' will select the corresponding mapping. 

Let us assume that we want to reach the same level of accuracy for all possible cases. 
As a consequence we observe that different processes under different sets of kinematical 
cuts will require values of NPTS and NRAND which may differ substantially. To illustrate 
this point we consider a particular example, e + e~ — > {i + fi~T + r~. Moreover we only 
include simple kinematical cuts on the invariant mass of the two fermion-antifermion 
pair, which can be coupled to a photon. The background to Z — Z production is therefore 
large, expecially when we allow for small cuts, M(/i + /i~), M(t + t~) > M with M much 
smaller than M z . To show our point we fix NPTS=8 and NRAND=6 and vary M . Table 
5 clearly illustrates as the numerical accuracy becomes better and better for growing M 
and fixed number of points. 



M (GeV) 


cr(fb) 


Rel. error(%) 


5 


10.29(8) 


0.77 


10 


6.61(3) 


0.43 


30 


3.130(2) 


0.07 


50 


2.226(1) 


0.03 



Table 5: The effect of varying the cut on fermion-antifermion pair in e + e — > t + t , 
inclusive of QED radiation. 

In this example we have used 0PEAK='N' for M = 5, 10 GeV and 0PEAK= 'Y' for M = 
30, 50 GeV. High precision for very loose cuts requires therefore a much higher number 
of points while it is considerably easier to reach 0.1% or better for more realistic cuts. 
It goes without saying that a more stringent cut on the uninteresting boundaries of the 
phase space implies high precision achievable with high computational speed. Finally we 
observe that for some of the processes with very loose cuts a repetition with an increasing 
sequence of high values of NPTS can eventually yield erratic results. 

ALS(REAL*8) 

For NC32 and non-leptonic NC64, Mix43 the value of a s (M w ) must be initialized. More- 
over for CC processes it is possible to include NQCD effects (Naive QCD). This is con- 
trolled by the variable OQCD. 

OQCD (CHARACTER* 1) 

Inclusion of NQCD both in the total W width and in the final state, ALS must be initial- 
ized. Thus 

r„ = C (i + Iv) ■ (120) 

0FL (CHARACTER*!) 
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Whenever a CC20 process is considered an option is available for introducing the gauge 
restoring terms according to the FL-scheme. The VF-width is always taken to be running. 
The e~ scattering angle cannot be set to zero even in the presence of gauge reparation 
due to the approximation of m e = 0. 

OFSR (CHARACTER* 1) 

In the present version of the program this flag has been set to its default, 0FSR='N' and 
it should not be changed. It is responsible for the inclusion of QED final state radiation 
for CC11 processes, a branch which deserves a more detailed theoretical analysis. 

ITCM ( INTEGER*4) 

For CC11 and CC20 several observables can be computed: 

• The cross-section 

• 1 unrenormalized moments of 

E, = (l - X -±±^-) VS. (121) 

• 2 for CC10. Unrenormalized moments of E^ as defined in eq . (|IT)|) . 

• 3 For CC10. Unrenormalized T n (cos where T n is the n-th Chebyshev polynomial 
and cos 9^ is defined in eq. (|H]) with i — 1. 

• 4 unrenormalized T n (cos9w-)- The W~ scattering angle is defined by 



cos 9 W 
E 

w 

Pi 



1 



1 r 

2 

1 



1 — x_ 



E 



(l + m 2 _ — m 2 ^ — (x + — X-) t l 



vm_ 
~EP~ 



(122) 



5 unrenormalized T n (cos9w+)- The W + scattering angle is defined by 



cos 9w+ 
E 



w+ 



E 



w+ 



(l+m+- ml) - (x + - a5_) (1 - t w ) 



0: 



w+ 



1 - 



w+ 



(123) 
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• 6 unrenormalized W virtuality defines as 



v w = Vv(m- +m + ) - 2 



w 



(124) 



s 



• 7 M(W + ) + M(W~) distribution 



• 8 M(W~) - M(W + ) distribution 



• 9 M{W + ) distribution 



• 10 c = cos#(/ ) distribution where 



t 1 = x + ei 



1 + c 



(125) 



x_ + x + — (ar_ — x + )c 



and theta is the angle between the incoming e and the outgoing I . 

while for all other classes only the cross-section will be available. Please consult the 
WW/eg report [|J for more detailed informations concerning this quantities. Here un- 
renormalized means not divided by the corresponding cross-section. 

ITCNM(INTEGER*4) 

If 1 < ITCM < 6 the order of the moments can be chosen. 
DIST(REAL*8) 

For ITCM=7,8,9 gives the value of M w + + M w -, M w - - M w + or M w + at which the 
differential cross section is required. For ITCM=10 DIST is 6l(e~,l~). 

Moreover for ITCM=0 there is the additional possibility of binning the M(W~) ± M{W + ) 
distributions by requiring the calculation of the cross-section with M 1 < M(W~) ± 
M(W+) < M 2 . 

OBIN (CHARACTER* 1) 

0BIN='N , ['PVM'] selects no binning or ± binning. For the latter we initialize 

ABP(M) ,BBP(M) (REAL*8) 
the limits M 1 ,M 2 . 

0C0UL (CHARACTER* 1) 

For CC processes the Coulomb correction factor can be included. 
I0S(INTEGER*4) 

Refers to the choice of the weak mixing angle, s 2 and of the SU(2) coupling constant g. 
The user is advised to adopt I0S=1 which is the choice of a large variety of programs and 
corresponds to (G M being the Fermi coupling constant) 
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V2G U M 2 ' 

f = ^ (126) 

w 



although a more sensible choice (connected with the Ward identities) is given by 



Ml 



w 

2 



M 2 ' 

z 



(127) 



ORAL ( CHARACTER* 1) 

If we adopt I0S=1 then the scale must be chosen for evaluating a. If 0RAL='F' then 
it is possible to enter any value for ALWI= 1/a while for 0RAL='R' the program will 
compute a(s) with a call to subroutine PSELF and HADR5. The two routines evaluate 
DER= DERL+DERH, DERL being the perturbative contribution to Ar from leptons and from 
the top quark while DERH evaluates the light hadron contribution using fits to the QED 
vacuum polarization from e + e~ data (courtesy of F. Jegerlehner 0). 

«W = T~^a~- (128) 
1 — Ar 



I0SF(INTEGER*4) 

Refers to the choice of the Structure Function for Initial State QED radiation. Born 
observables require I0SF=0, for the corrected ones the user is advised to have I0SF=1. 
Indeed this choice is the default adopted by the WW/eg working group and corresponds 
to having D(x, s) of eq.( |97|) entirely in terms of (3 ||, 

D r (x, s) = ~ (1 - *Y m (l ~x 2 )+0 ((3) . (129) 
Instead I0SF=2 corresponds to the so-called //-scheme [|J 




Dr(x,8) = -\^-xY m 

01 Q 

V = 2- log—, (130) 

which respects gauge invariance. Finally I0SF=3 gives a mixed treatment as described 
in [0. There is no way in which one can present any realistic estimate of the theoreti- 
cal error for 4-fermion production with the actual level of knowledge that we have. The 
only pragmatic alternative is to examine the variation induced in our results by changing 
RS and SF. The first option roughly simulates (probably underestimate) the uncertainty 



33 



connected with our ignorance of the 0(a) electroweak corrections while the second one 
reflects the uncertainty related to our ignorance about the complete 0(ax const) QED 
corrections which, at the moment, only control correctly the leading logarithmic parts. 
This we have done by allowing I0S=1, 2 and I0SF=1, 2, 3 in the calculation of the repre- 
sentative process e + e~ — > n~V^, v t t + . Results for the cross section at y/s = 190 GeV and 
for canonical cuts are show in table 6. 



I0S/I0SF 


1 


2 


3 


1 
2 


0.19411(2) 
0.19394(2) 


0.19395(2) 
0.19378(2) 


0.19464(2) 
0.19447(2) 



Table 6: a (nb) for e + e — > /i z/ M , u t t + at 190 GeV with E^ T > 1 GeV and 
10° < e^ T < 170°. Moreover 9^ T > 5°. 

Thus naively one could fix the central value at I0S=2 and I0SF=2 which preserve gauge 
invariance and assign a quite asymmetric theoretical error of 

a = 0.19378 ± 0.00002(num.)+°;™(theor.). (131) 

Moreover we have also show in tables 7 — 8 the effect of varying I0S/I0SF in the process 
e + e~ — > fi'u^ud both for a and for < >. In this case we have included the Coulomb 
correction factor, 0C0UL='Y', and the NQCD corrections, 0QCD='Y'. Canonical cuts are 
applied. The largest uncertainty for < E 1 > is of 1.9, 3.2, 9.9, 20.5 MeV for E cm = 
161, 175, 190, 205 GeV respectively. 



E cm (GeV)/IOS-IOSF 


1-1 


1-2 


1-3 


2-1 


2-2 


2-3 


161 


0.4688 


0.4673 


0.4674 


0.4685 


0.4669 


0.4670 


175 


1.1250 


1.1219 


1.1221 


1.1251 


1.1220 


1.1222 


190 


2.1579 


2.1484 


2.1489 


2.1583 


2.1488 


2.1493 


205 


3.2317 


3.2119 


3.2129 


3.2324 


3.2126 


3.2135 



Table 7: < £ 7 > (GeV) for e+e' -> yTv^ud with E^ > 1 GeV, E u4 > 3 GeV, 
M ud > 5 GeV, 10° < 0„ < 170°, 9^ u , 9^ d > 5°. 

0CUTS(CHARACTER*2) It is the main decision branch for kinematical cuts. 
There are few internal options. 

1. 0CUTS='EX'. Fully extrapolated set-up, use it only for processes free of singularities, 
e.g. never for CC20. 

2. 0CUTS='CC. The Canonical Cuts used within the WW/eg Working Group. In 
this case cuts are automatically selected. Canonical Cuts are also available for 
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E cm (GeV)/IOS-IOSF 


1-1 


1-2 


1-3 


2-1 


2-2 


2-3 


161 


0.13206 


0.13204 


0.13250 


0.13201 


0.13198 


0.13244 


175 


0.49207 


0.49186 


0.49358 


0.49177 


0.49156 


0.49329 


190 


0.60240 


0.60192 


0.60404 


0.60188 


0.60139 


0.60352 


205 


0.62828 


0.62754 


0.62977 


0.62764 


0.62691 


0.62913 



Table 8: a (nb) for e + e — > fi v^ud. The cuts are the same as in table reftabl. 

NC21,NC25 and NC50 processes. For completeness we make a list of the chosen 
canonical cuts ||. Let I be any final state charged lepton and q any final state 
quark. Then 

• Ei > 1 GeV, E q > 3 GeV. 

• Mfaqj^MfaqjlMfaqj) > 5 GeV. 

• 10° <9i< 170°. 

• e^q^^ik,^) >5°. 

3. 0CUTS='SC. No cuts at all but whenever a // pair appears in the final state a cut 
is applied to prevent singularity due to 7* — > //. Again not to be used for CC20 
etc. Actually it can ONLY be used for CC11, NC19(NC21), NC24(NC25), NC50 
processes. The 'security' cut for CC11, NC19 and NC24 can be initialized through 
the variable STHG. For Higgs physics, NC21, NC25 and NC50 processes, the 'simple' 
cuts can be initialized through the variables ZCUT , BCUT. 

4. 0CUTS='FC. You adopt your own cuts. Several variables must be initialized. The 
important thing to pay attention to is the 'ordering' of the final states. One should 
have in mind the classification given in appendix, therefore if IPR=18 is selected 
one must attribute cuts according to the fact that 1 = u, 2 = u, 3 = u e , 4 = V e . The 
internal coding will attribute the following labels: for CC(Mix) processes e + e~ — > 
d(l)u(2)u\Z)d' (4) , where u = fi,u,c and d = l,d,s,b. For NC processes e + e _ — > 
/(l)/(2)/'(3)/ (4). Boundaries for the scattering angles or final state energies are 
stored as 1 — 4. Boundaries for final state angles 0^, i, i = 1, 4 or final state invariant 
masses My, i, % = 1, 4 are stored as a vector of dimension 6, with entries 1 — 2/1 — 
3/1 -4/2-3/2-4/3-4. 

STHG(REAL*8) 

Selects M(ff) > STHG. 

ZCUT, BCUT (REAL*8) 

In e+e" -> Jfbb selects M z - ZCUT < M(ff) < M z + ZCUT and M(bb) > BCUT. 
AIMM(I) (REAL*8) 
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Min cuts in GeV on FS invariant masses 1 — 6. 
BIMM(I) (REAL*8) 

Max cuts in GeV on FS invariant masses 1—6. 
AEM(I) (REAL*8) 

Lower cuts on FS energies (GeV) 1—4. 
ASAM(I) (REAL*8) 

Min cuts on the cosine of the scattering angles 1 — 4. 
BSAM(I) (REAL*8) 

Max cuts on the cosine of the scattering angles 1—4. 
AFSAM(I) (REAL*8) 

Min cuts on the cosine of the FS angles 1 — 6. 
BFSAM(I) (REAL*8) 

Max cuts on the cosine of the FS angles 1 — 6. 



Whenever the Higgs signal has to be investigated few additional flags must be initialized. 
They are: 



HMI (REAL*8) 

The Higgs boson mass in GeV. 



ALSM(REAL*8) 

The value of a s (M w ). 



0QCD( CHARACTER* 1) 

The option for final state 'naive' QCD corrections. 



0GLU( CHARACTER* 1) 

The inclusion of the H — > gg channel in the total Higgs width. The most complete 
treatment will therefore evolve a s to the scale \i = M H , evaluate the running b, c-quark 
masses and compute 



G G M B 
An 



m 2 b (M H )+m 2 c (M H ) 



a q tot, 
1 + 5.67— + 42.74 — 

n \n 



+ m + r 



99' 



99 



^fi4f 1+17 . 918e7 «. 

d6/T n l \ n 



(132) 



Again the 'ordering' of final states is extremely important when one chooses cuts. Please 
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notice that for internal reasons WTO attributes the following labels to IPR=30-33: 1 = 

6,2 = 6,3 = ui/l-,4 = Vi/l + but 1 = fi~, 2 = 3 = b, 4 = 6. 

2.10 Test run output 

The typical calculations that can be performed with the program are illustrated in three 
different examples. 

Sample 1. Evaluation of the total cross section for the process e + e~ — > fx~, u, d inclu- 
sive of QED initial state radiation. 

Sample 2. Evaluation of the total cross section for the process e + e~ — > e + e~dd without 
QED initial state radiation. 

Sample 3. Evaluation of the total cross section for the process e + e~ — > /i + /!~z7 M z/ M with 
QED initial state radiation. 

In sample 1 - 3 we choose E cm = 190 GeV, M w = 80.23 GeV, M z = 91.1888 GeV and 
T z = 2.4974 GeV. The various options are as follows. 
For sample 1 we have: 

0PEAK= , Y ) 

The double-resonating mapping is requested. 
ITC=0 

The total cross section is computed. 
DC0UL='N' 

The Coulombic correction factor is not included. 
IOS = 1, IOSF = 1 

The default choices for the renormalization scheme and the structure functions. 

0RAL= ' F ' , ALWI = 128.07D0 

We have fixed a to be l/a = 128.07/JO. 

OCUT = 'CC 

Canonical Cuts are applied. 
IPR = 2 

The internal process number is selected. 
For sample 2 we have: 
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OPEAK='N' 

The double-resonating mapping is NOT requested. 
0C0UL='N' 

The Coulombic correction factor does not apply for this process. 
IOS = 1, IOSF = 

The default choice for the renormalization scheme and NO QED radiation. 

0RAL= ' F ' , ALWI = 128.07D0 

We have fixed a to be l/a = 128.07D0. 

OCUT = 'CC 

Canonical Cuts are applied. 
IPR = 25 

The internal process number is selected. 

For sample 3 we have: 
0PEAK='Y' 

The double-resonating mapping is requested. 
0C0UL='N' 

The Coulombic correction factor is NOT not requested. 
IOS = 1, IOSF = 1 

The default choices for the renormalization scheme and the structure functions. 

0RAL= ' F ' , ALWI = 128.07D0 

We have fixed a to be l/a = 128.07D0. 

OCUT = 'CC 

Canonical Cuts are applied. 
IPR = 20 

The internal process number is selected. 
2.11 Program implementation 

In this section a short description of each routine implemented in WTO is given. 

BLOCK DATA INIT 

IMPLICIT REAL* 8 (A-H.O-Z) 
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CHARACTER* 1,0PRT 

* 

COMMON/PRT/OPRT 
COMMON/PARAM/EPS , QDELTA 

COMMON/FMASS/EM , RMM , TM , TQM , UQM , DQM , CQM , SQM , BQM 
COMMON/BPAR/WM , ZM , ZG , GF , PI , CFCT , FCNT , GE , ALPHAI , ALWI 

* 

* OFFICIAL SET OF DATA INPUT 

* 

DATA WM/80 . 23D0/ , CFCT/O . 38937966D9/ , FCNT/81 . DO/ , 

# GF/8 . 2476172696D-6/ , ZM/91 . 1888D0/ , 

# EM/O . 51099906D-3/ ,PI/3 . 141592653589793238462643D0/ , 

# GE/O . 5772156649D0/ , ZG/2 . 4974D0/ ,TM/1 . 7771D0/ , 

# BQM/4 . 7D0/ , CQM/1 . 55D0/ , ALPHAI/137 . 0359895D0/ 
DATA EPS/1. D-37/,QDELTA/0. DO/ 

DATA RMM/O . 10565839D0/ ,UQM/0 . 041DO/ ,DQM/0 . 041D0/ , SQM/O . 15D0/ 
DATA TQM/175. DO/ 
DATA OPRT/'Y'/ 

* 

END 



Masses and numerical constants are initialized. As it can be seen the top quark mass is 
also initialized. At the moment this is only due to the following facts. Whenever a is 
evoluted at the scale s then the top quark enters the photon two point function. Morever 
in NC21, i.e. e + e~ — > bbv e V e) there is a multiperipheral diagram whith an internal top 
line. 



SUBROUTINE TOW (RS , I ACM , AIMM , BIMM , A EM , ISAAM , IS ABM , ASAM , BSAM , 

# AFSAM , BFSAM , OTYPEM , CHUM , CHUPM , CHDM , CHDPM , CHFM , 

# CHFPM , TIFM , TIFPM , FCUM , FCDM , ITCM , IPM , IRM , 

# ITCNM.OXCM) 
IMPLICIT REAL* 8 (A-H.O-Z) 

CHARACTER* 1 , OXCM , OQCD , OGLU , OFL , OCOUL , OFSR , OTRANS , OPRT 
CHARACTER*2 , OFS 
CHARACTER* 5 RLABS (3) 

CHARACTER* 7 CLABS1 (4) ,CLABS2(4) ,CLABS3(4) ,CLABS4(4) , 



# CLABS5 (4) , CLABS6 (4) , CLABS7 (4) , CLABS8 (4) , 

# CLABS9 (4) ,CLABS10(4) ,CLABS11(4) ,CLABS12(4) , 

# CLABS13 (4) ,CLABS14(4) ,CLABS15(4) ,CLABS16(4) , 

# CLABS17 (4) ,CLABS18(4) ,CLABS19(4) ,CLABS20(4) , 

# CLABS21 (4) ,CLABS22(4) ,CLABS23(4) ,CLABS24(4) , 

# CLABS25 (4) , CLABS26 (4) , CLABS27 (4) , CLABS28 (4) , 

# CLABS29 (4) , CLABSO (4) 



CHARACTER*4 , OTYPEM , OTYPE 

PARAMETER (N0UT=99 , NDIMMX=9) 

COMMON/MP/ZRM 
COMMON/FS/OFS 
COMMON/GG/OGLU 
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COMMON/FLS/OFL 

COMMON/QCD/ALS 

COMMON/NPR/IPR 

COMMON/HIGGS/HM 

COMMON/PRT/OPRT 

COMMON/FSR/OFSR 

COMMON/OPA/DELC 

COMMON/DIS/DIST 

COMMON/KOUNT/IK 

COMMON/SP/PSGC4) 

COMMON/ISTRF/ISF 

COMMON/AQCD/OQCD 

COMMON/BME/BFACT 

COMMON/COUL/OCOUL 

C0MM0N/CHI/HCH(36) 

COMMON/TC/ITC,ITCN 

COMMON/SHEL/OTRANS 

COMMON/IPT/IFZ(44) 

COMMON/ICUTS/IAC (4) 

C0MM0N/RM/RBM2 , RCM2 

COMMON/TOPT/IOS , IDSF 

COMMON/ISA/ISAA, ISAB 

COMMON/OCHANNEL/OTYPE 

COMMON/NCC/ CHF2 , CHFP2 , CONC ( 10) 

COMMON/CCLR/VUPL , VUPR , VDPL , VDPR 

COMMON/NCLR/VEL , VER , VELR , VFL , VFR , VFPL , VFPR 

COMMON/CCHANNEL/CHU, CHUP , CHD , CHDP , FCUC , FCDC 

COMMON/NCHANNEL/CHF , CHFP , TIF , TIFP , FCUN , FCDN 

COMMON/FMASS/EM , RMM , TM , TQM , UQM , DQM , CQM , SQM , BQM 

COMMON/HAP AR/RHM , RHM2 , RHG , RHMG , SHG , SHGS , OPSHGS 

C0MM0N/EE/QCH,QCH2,VQR,VQL,HBE(24) ,HB0(24) ,HMP(24) 

COMMON/BPAR/WM , ZM , ZG , GF , PI , CFCT , FCNT , GE , ALPHA I , ALWI 

COMMON/CPAR/ ALPHA , HBET , HBETI , OMHB , EDB , DOGL , G8 , TFACT , PIH , ALW , 



# ETA, FETA, BETA, G2 

COMMON/ ACCHANNEL/ DMCHU , OPCHU , OMCHUP , OPCHUP , OMCHDP , OPCHDP , 

# OMCHD , OPCHD , HCHUP , HCHU , HCHDP , HCHD 
COMMON/APAR/ARS , S , RWM , RWM2 , RWG , RWMG , SWG , SWGS , OPSWGS , STH2 , CTH2 , 

# HSTH2 , TSTH2 , SCTH2 , ASTH2 , TTH2 , RZM , RZM2 , RZG , RZMG , SZG , 

# SZGS , OPSZGS , STH4 , CTH4 , VE , VF , VFP , RBQM2 , RSZW , RSZW2 , 

# SOW.SOZ 

COMMON/ SUBREG/DSM , USM , DSP , USP , RL ( 6 ) , RR ( 6 ) , SRL ( 6) , SDSM , SDSP , VVL1 , 

# VVL2,VVL3,UL,0MUL,SUML 

COMMON/CUTS/AIM (6) ,BIM(6) ,AE(4) ,ASA(4) ,BSA(4) ,AFSA(6) ,BFSA(6) , 

# OMBSA (4) , OPBSA (4) , TEQ , RAE (4) , OMASA (4) , OPASA (4) , 

# SG12 , CG12 , SG13 , CG13 , SG14 , CG14 , SG23 , CG23 , SG24 , 

# CG24 , SG34 , CG34 , SCT120 , SCT130 , SCT140 , SCT230 , 

# SCT240,SCT340,SGAM(4) ,CGAM(4) 



DIMENSION VK(NDIMMX) 

DIMENSION IACM(4) ,AIMM(6) ,BIMM(6) ,AEM(4) ,ASAM(4) ,BSAM(4) , 

# AFSAM(6) ,BFSAM(6) ,CUTS1(3,4) ,CUTS2(4,4) ,CUTS3(4,4) , 

# CUTS4 (4,4), CUTS 5 (4,4) 
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EXTERNAL D01GCF , REGION , XSC , XSN , XSM43 , XSNH19 , XSNH24 , XSNH49 , XSN32 , 
# XSN64 , S14AAF , S21BAF , X04CBF , X04ABF , S09ABF 

This subroutine will decide, according to the input parameters, which process to consider, 
e.g. which class and family, and which observable to compute. For each diagram all the 
couplings are computed, i.e. given a current 

7"(a + & 7 5 ), (133) 

the chiral couplings constants 

9L = \(" + b), 9 R = \ip,-b), (134) 

are stored and all the relevant products Hig 1 n^g- 7 are computed once for all. 

All the basic parameters which control the boundaries of the phase space are initialized 
in TOW. All the I/O is controlled by this subroutine after a call to the subroutine D01GCF 
which performs the integration. The following functions give the differential cross sections 
according to the family to which the process belong. All of them essentially set the 
ingredients for the phase space and the requested mappings. Once the vector X(NDIM) is 
initialized, with < Xi < 1, i — 1, . . . , NDIM, all independent invariants u, v, m 2 _, . . . , t l are 
constructed and from them all the combinations Xij, i < j = 1, . . . , 6. Next all non-zero 
helicity amplitudes are computed and squared. Within these subroutines the operations 
among complex numbers are always executed in terms of real and imaginary parts. 

REAL* 8 FUNCTION XSC (NDIM, X) 
IMPLICIT REAL* 8 (A-H.O-Z) 
CHARACTER* 1 , OCOUL , OPEAK , OQCD , OFL 
CHARACTER* 10,0TYPE 
DIMENSION X(NDIM) 

With XSC all the CC3-CC11-CC20 processes are computed. As explained in the long 
write-up this function will allow different choices for the observables to be computed, 
including the total cross section and the moments of several distributions. Moreover the 
variable 0FL='Y[N]' controls the inclusion of the imaginary part of the VWW, V = 7, Z 
one loop vertex in CC20. 

REAL* 8 FUNCTION XSNH24(NDIM,X) 
IMPLICIT REAL* 8 (A-H.O-Z) 
CHARACTER* 1, OQCD 
COMMON/HIGGS/HM 
COMMON/AQCD/OQCD 
DIMENSION X(NDIM) 
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With this function the NC25 processes are compute by means of calls to functions XSH24 
and XSN. If 0QCD='Y' then a s (M H ) is computed through a call to function RALPHAS and 
corrections are applied. 

REAL* 8 FUNCTION XSNH19(NDIM,X) 
IMPLICIT REAL* 8 (A-H.O-Z) 
CHARACTER* 1,0QCD 
COMMON/AQCD/OQCD 
DIMENSION X(NDIM) 



The same as for function XSNH24 but for NC21 processes. 



REAL* 8 FUNCTION XSNH49(NDIM,X) 
IMPLICIT REAL* 8 (A-H.O-Z) 
CHARACTER* l.OQCD 

* 

DIMENSION X(NDIM) 



The same as for function XSNH24 but for NC50 processes. 



REAL* 8 FUNCTION XSN32(NDIM,X) 
IMPLICIT REAL* 8 (A-H.O-Z) 
DIMENSION X(NDIM) 

This function will compute all NC32 processes by first calling function XSN, i.e. the NC24 
sub-class (no gluons), and then by calling function XSNG which will add the remaining 
8 diagrams. This is made possible by the fact that there is no interference between NC 
diagrams and QCD diagrams. For historical reasons the corresponding Born cross sections 
(I0SF=0) are computed without gluon exchange diagrams for non-leptonic processes. 

REAL* 8 FUNCTION XSM43 (NDIM , X) 
IMPLICIT REAL* 8 (A-H.O-Z) 
CHARACTER*2 , OFS 
COMMON/FS/OFS 
DIMENSION X(NDIM) 

This function will compute both leptonic and non-leptonic Mix43 processes by first calling 
function XS35, i.e. the 35 diagrams without gluons but with both CC and NC contribu- 
tions. Next function XS35G is called which includes the QCD diagrams and the interfer- 
ence CCcg>QCD. For historical reasons the corresponding Born cross sections (I0SF=0) 
are computed without gluon exchange diagrams for non-leptonic processes. 
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REAL* 8 FUNCTION XSN(NDIM,X) 
IMPLICIT REAL*8 (A-H.O-Z) 
CHARACTER* l.OPEAK 
CHARACTER* lO.OTYPE 
COMMON/PS/OPEAK 
COMMON/OCHANNEL/ OTYPE 
DIMENSION X(NDIM) 

With this function all available NC processes are computed but without Higgs exchange. 
Thus it includes the NC19, NC24 and NC48 families. Only the total cross section is 
available. 

REAL* 8 FUNCTION XSH24(NDIM,X) 
IMPLICIT REAL* 8 (A-H.O-Z) 
DIMENSION X(NDIM) 

This function will include the Higgs boson exchange diagram to the NC24 processes. It 
should be stressed that in the limit of massless fermions there is no interference between 
NC24 and Higgs so that one can compute signal and background separately. If 0QCD='Y' 
then a s (M w ) is evoluted to the a s (M H ) and the Higgs boson signal is multiplied by 

5 QCD = 1 + 5.67^ + 42.74 * , 

a s = a s (M H ). (135) 



REAL* 8 FUNCTION XSH19(NDIM,X) 
IMPLICIT REAL* 8 (A-H.O-Z) 
DIMENSION X(NDIM) 

With this function the Higgs boson exchange is added to the NC19 family thus obtaining 
NC21. The two additional diagrams again do not interfere with NC19 in the massless 
limit. 

REAL* 8 FUNCTION XSNH49(NDIM,X) 
IMPLICIT REAL*8 (A-H.O-Z) 
CHARACTER* l.OQCD 

* 

COMMON/QCD/ALS 

COMMON/HIGGS/HM 

COMMON/AQCD/OQCD 

DIMENSION X(NDIM) 
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With this function the Higgs boson exchange is added to the NC48 family thus obtaining 
NC50. The two additional diagrams again do not interfere with NC48 in the massless 
limit. 

REAL* 8 FUNCTION XSNG(NDIM,X) 
IMPLICIT REAL* 8 (A-H.O-Z) 
CHARACTER* 10,0TYPE 
COMMON/OCHANNEL/DTYPE 
DIMENSION X(NDIM) 

The gluon exchange diagrams for non-leptonic NC processes are included. As explained 
in the previous sections we use a fixed a s , i.e. a s (M w ). 

REAL* 8 FUNCTION XS35(NDIM,X) 
IMPLICIT REAL* 8 (A-H.O-Z) 
CHARACTER* 1 , OCOUL , OPEAK 
CHARACTER*2 , OFS 
CHARACTER* 10,0TYPE 
COMMON/FS/OFS 
COMMON/PS/OPEAK 
C0MM0N/C0UL/0C0UL 
COMMON/OCHANNEL/OTYPE 
DIMENSION X(NDIM) 

This function computes the diagrams for CC®NC leptonic processes and for the non- 
leptonic ones where gluon exchange will be included separately. 

REAL* 8 FUNCTION XS35G(NDIM,X) 
IMPLICIT REAL* 8 (A-H.O-Z) 
CHARACTER* 1 , OCOUL , OPEAK 
CHARACTER*2 , OFS 
CHARACTER* 10,0TYPE 
COMMON/FS/OFS 
COMMON/PS/OPEAK 
C0MM0N/C0UL/0C0UL 
COMMON/OCHANNEL/OTYPE 
DIMENSION X(NDIM) 



Whenever CC®NC non-leptonic processes are to be computed this function will return 
QCD®(QCD+CC). 

REAL* 8 FUNCTION XSN64(NDIM,X) 
IMPLICIT REAL*8 (A-H.O-Z) 
CHARACTER* 1, OPEAK 
CHARACTER*2 , OFS 
COMMON/FS/OFS 
COMMON/PS/OPEAK 
DIMENSION X(NDIM) 
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Finally this function allows to compute processes with identical particles in the final state, 
with the exclusion of and of v e . For the non-leptonic case QCD is included and the 
proper color factor relative to each diagram is automatically computed. 

SUBROUTINE REGION (N ,X , J , A ,B) 
IMPLICIT REAL*8(A-H,0-Z) 
DIMENSION X(N) 

This subroutine is requested by D01GCF in order to set the limits of integrations. In 
WTO its content is trivial since all variables are in the range [0, 1]. 

SUBROUTINE CORRQCD (SCAL , ALS , BQM2 , COM2 , RBM2 , RCM2) 
IMPLICIT REAL*8(A-H,0-Z) 

COMMON/BPAR/WM , ZM , ZG , GF , EM , PI , CFCT , FCNT , GE , TM , BOM , COM , ALWI 
DATA Z3/1.20206D0/,NF/5/ 

Computes the QCD quark running masses. Given ALS, the value of a s (M w ), and the pole 
b, c-quark masses this routine evaluates their running at a scale /i = SCAL. 

REAL* 8 FUNCTION RALPHAS (RSO ,RS , ALS ,NF) 
IMPLICIT REAL*8(A-H,0-Z) 

COMMON/BPAR/WM , ZM , ZG , GF , EM , PI , CFCT , FCNT , GE , TM , BOM , COM , ALWI 

Given ALS, the value of a s at a scale RSO, this routine returns a s at the scale [i = RS with 
NF active flavors. 

REAL* 8 FUNCTION QCDLAM (NF , ALS , RS , XI , X2 , XACC) 
IMPLICIT REAL*8(A-H,0-Z) 
PARAMETER ( JMAX=50 , N0UT=21) 

REAL* 8 FUNCTION QCDSCALE (NF , ALS , RS , X) 
IMPLICIT REAL*8(A-H,0-Z) 



Compute A n -f in the MS'-scheme from a s fixed at the scale RSO. 



SUBROUTINE HADR5 (E , DER , EDER) 



* * 

* SUBROUTINE FOR THE EVALUATION OF THE LIGHT HADRON * 

* CONTRIBUTIONS TO DELTA_R * 

* * 

* F. JEGERLEHNER, PAUL SCHERRER INSTITUTE, CH-5232 VILLIGEN * 

* * 
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* REFERENCE: F. JEGERLEHNER , Z. PHYS . C32 (1986) 195 * 

* H. BURKHARDT ET AL., Z. PHYS. C42 (1989) 497 * 

* S. EIDELMAN , F. JEGERLEHNER, Z. PHYS. C (1995) * 

* * 

* VERSION: 24/02/1995 
* 

* THIS PROGRAM DOES NOT YET KNOW HOW TO COMPUTE DELTA R AND DELTA G FOR 

* ENERGIES IN THE RANGES |E|>1TEV AND 2M_PI < E < 40(13) GEV !!!!!!!!! 

IMPLICIT REAL*8(A-H,0-Z) 



SUBROUTINE PSELF(P2X,PGGF) 
IMPLICIT REAL*8(A-H,0-Z) 

COMMON/PARAM/EPS , QDELTA 

COMMON/FMASS/EM , RMM , TM , TQM , UQM , DOM , COM , SQM , BOM 
COMMON/BPAR/WM , ZM , ZG , GF , PI , CFCT , FCNT , GE , ALPHAI , ALWI 



PSELF computes the perturbative contribution to the running of a(s). 



SUBROUTINE RBFFO (RM12 , RM22 , BO , B 1 , B2 1 ) 
IMPLICIT REAL*8(A-H,0-Z) 

COMMON/PARAM/EPS , QDELTA 

COMMON/BPAR/WM , ZM , ZG , GF , PI , CFCT , FCNT , GE , ALPHAI , ALWI 

DIMENSION ARG(2) ,CLN(2) ,FR(3) 

SUBROUTINE RBFF (P2 , RM12 , RM22 , RBO , RBI , RB2 1 ) 
IMPLICIT REAL*8(A-H,0-Z) 

COMMON/PARAM/EPS , QDELTA 

COMMON/BPAR/WM , ZM , ZG , GF , PI , CFCT , FCNT , GE , ALPHAI , ALWI 

DIMENSION CMP2(2) 
DIMENSION CLNMP2(2) 
DIMENSION GFPR(3) ,GFMR(3) 

SUBROUTINE ROOTS (P2 , RM12 , RM22 , RPR , RPI , RMR , RMI , OMRPR , OMRMR) 
IMPLICIT REAL*8(A-H,0-Z) 

COMMON/PARAM/EPS , QDELTA 

COMMON/BPAR/WM , ZM , ZG , GF , PI , CFCT , FCNT , GE , ALPHAI , ALWI 

SUBROUTINE CQLNX(ARG.RES) 
IMPLICIT REAL*8(A-H,0-Z) 

COMMON/BPAR/WM , ZM , ZG , GF , PI , CFCT , FCNT , GE , ALPHAI , ALWI 
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DIMENSION ARG(2) ,AARG(2) ,RES(2) 

SUBROUTINE CQLNOMX (ARG , OMARG , RES) 
IMPLICIT REAL*8(A-H,0-Z) 

DIMENSION ARG (2) , OMARG (2) ,RES(2) , ARES (2) ,CT(10) ,SN(10) 

SUBROUTINE RCG (N , ZR , ZI , OMZR , GFR) 
IMPLICIT REAL*8(A-H,0-Z) 

COMMON/PARAM/EPS , QDELTA 

COMMON/BPAR/WM , ZM , ZG , GF , PI , CFCT , FCNT , GE , ALPHA I , ALWI 
DIMENSION GFR(N) 

DIMENSION A(4) , RA (16) , Z (2) , 0MZ(2) , OZ (2) ,CLN0MZ(2) ,CLN0Z(2) 
DIMENSION CA(2,10) ,AUX(2) ,ZP(2) ,CT(16) ,SN(16) 



These routines are required by PSELF and compute the relevant perturbative contributions 
to the photon two point function. 



3 Test run output 
3.1 Sample 1 

This run is with: 

NPTS = 8 

NRAND = 6 

l/alpha_QED(s) = 0.12807E+03 Coupling constants are (alpha-G_F) : 



gv_e 
ga_e 
g-wf 

g_ZWW 

G-L factor = 

E_cm (GeV) = 
beta 

M_W (GeV) = 

G_W (GeV) = 



alpha RS 

0(alpha~2) beta SF 
CS not included 
FSR not included 
NQCD not included 



-0.1409725455532E-01 
-0.1857939837526E+00 
0.2304098503535E+00 
-0.5714792058132E+00 

. 104114673716395E+01 



-0.1409725455532E-01 
-0.1857939837526E+00 
0.2304098503535E+00 
-0.5714792058132E+00 



0.19000E+03 
0.11453E+00 sin~2 
0.80230E+02 M_Z (GeV) = 
0.20337E+01 G_Z (GeV) = 



0.23103E+00 
0.91189E+02 
0.24974E+01 



Process is e+ e- -> mu- 



bnu_mu 



bd 
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with the following cuts 



E(GeV), SA(deg) cuts 



mu- bnu_mu u 

E_th 1.0000 0.0000 3.0000 3, 

C_max 10.0000 0.0000 0.0000 0. 

C_min 170.0000 180.0000 180.0000 180. 

lower IM cuts (GeV) 



bd 
0000 
0000 
0000 



mu- 
bnu_mu 
u 
bd 



bnu_mu u bd 

0.0000 0.0000 0.0000 
0.0000 0.0000 
5.0000 



upper IM cuts (GeV) 



mu- 
bnu_mu 
u 
bd 



bnu_mu u 
190.0000 190.0000 190. 

190.0000 190. 

190 



bd 

0000 
0000 
0000 



min FS angle (deg) cuts 



mu- 
bnu_mu 
u 
bd 



mu- bnu_mu u 

180.0000 180.0000 180 
180.0000 180 
180 



bd 
0000 
0000 
0000 



max FS angle (deg) cuts 



mu- 



mu- 
bnu_mu 
u 



bnu_mu 
0.0000 



u 

5.0000 
0.0000 



bd 
0000 
0000 
0000 



cell-diagrams : charges -1.0000 0.0000 

Cross-Section On exit IFAIL = 

CPU time 4 h 15 min 5 sec 

sec per call = 0.249E-02 

sigma = 0.5919248E+00 +- . 4040678E-04 

Rel. error of 0.007 '/, 



0.6667 



0.3333 
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3.2 Sample 2 



This run is with: 



NPTS 
NRAND 



10 
5 



l/alpha_QED(s) = 0.12807E+03 Coupling constants are (alpha-G_F) 



gv_e 
ga_e 
g_wf 

g_ZWW 

G-L factor 



-0.1409725455532E-01 
-0.1857939837526E+00 
0.2304098503535E+00 
-0.5714792058132E+00 



-0.1409725455532E-01 
-0.1857939837526E+00 
0.2304098503535E+00 
-0.5714792058132E+00 



. 104114673716395E+01 



E_cm (GeV) = 
beta 

M_W (GeV) = 

G_W (GeV) = 



0.19000E+03 
0.11453E+00 sin~2 
0.80230E+02 M_Z (GeV) 
0.20337E+01 G_Z (GeV) 



0.23103E+00 
0.91189E+02 
0.24974E+01 



alpha RS 

No QED Radiation 



Process is e+ e- -> d 
with the following cuts 



bd 



C_min 



e+ 



E(GeV) , SA(deg) cuts 

d bd e- e+ 

E_th 3.0000 3.0000 1.0000 1.0000 

C_max 0.0000 0.0000 10.0000 10.0000 



180.0000 180.0000 170.0000 170.0000 



lower IM cuts (GeV) 

d bd 

d 5.0000 
bd 
e- 
e+ 

upper IM cuts (GeV) 

d bd 

d 190.0000 
bd 
e- 
e+ 

min FS angle (deg) cuts 

d bd 

d 180.0000 
bd 
e- 
e+ 





e- 




e+ 


0. 


.0000 


0. 


0000 


0. 


.0000 


0. 


0000 






0. 


0000 




e- 




e+ 


190. 


.0000 


190. 


.0000 


190. 


,0000 


190. 


.0000 






190. 


.0000 




e- 




e+ 


180. 


,0000 


180. 


0000 


180. 


.0000 


180, 


.0000 






180. 


.0000 
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max FS angle (deg) cuts 

d bd 

d 0.0000 
bd 
e- 
e+ 

nc48-diagrams : charges 
isospin 



e- 

5.0000 
5.0000 



-0.3333 
-0.5000 



e+ 

5.0000 
5.0000 
5.0000 



-1.0000 
-0.5000 



Cross-Section On exit IFAIL = 

CPU time 9 h 9 min 3 sec 

sec per call = 0.111E-02 

sigma = 0.4304695E-01 +- . 5671371E-04 

Rel. error of 0. 132 7, 



3.3 Sample 3 



This run is with: 



NPTS 
NRAND 



l/alpha_QED(s) = 0.12807E+03 Coupling constants are (alpha-G_F) 



gv_e 
ga_e 
g-wf 
g_zww 



-0.1409725455532E-01 
-0.1857939837526E+00 
0.2304098503535E+00 
-0.5714792058132E+00 



-0.1409725455532E-01 
-0 . 1857939837526E+00 
0.2304098503535E+00 
-0.5714792058132E+00 



G-L factor 



. 104114673716395E+01 



E_cm 
beta 
M_W 
G_W 



(GeV) 

(GeV) 
(GeV) 



0.19000E+03 
0.11453E+00 sin~2 
0.80230E+02 M_Z (GeV) 
0.20337E+01 G_Z (GeV) 



0.23103E+00 
0.91189E+02 
0.24974E+01 



alpha RS 

0(alpha~2) beta SF 
CS not included 



Process is e+ e- -> mu- 
with the following cuts 



bnu_mu nu_mu mu+ 



E(GeV) , SA(deg) cuts 

mu- bnu_mu nu_mu mu+ 

E_th 1.0000 0.0000 0.0000 1.0000 
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C_max 
C_min 



10.0000 
170.0000 



0.0000 
180.0000 



0.0000 
180.0000 



10.0000 
170.0000 



lower IM cuts (GeV) 

mu- bnu_mu 
mu- 0.0000 
bnu_mu 
nu_mu 
mu+ 



nu_mu mu+ 
0.0000 0.0000 
0.0000 0.0000 
0.0000 



upper IM cuts (GeV) 

mu- bnu_mu nu_mu 
mu- 190.0000 190.0000 

bnu_mu 190.0000 
nu_mu 
mu+ 

min FS angle (deg) cuts 

mu- bnu_mu nu_mu 
mu- 180.0000 180.0000 

bnu_mu 180 . 0000 

nu_mu 
mu+ 

max FS angle (deg) cuts 

mu- bnu_mu nu_mu 
mu- 0.0000 0.0000 

bnu_mu . 0000 

nu_mu 
mu+ 

Cross-Section On exit IFAIL = 

CPU time 7 h 59 min 2 sec 

sec per call = 0.139E-02 

sigma = 0.2046758E+00 +- . 1704597E-03 

Rel. error of 0.083 7„ 



mu+ 
190.0000 
190.0000 
190.0000 



mu+ 
180.0000 
180.0000 
180.0000 



mu+ 
5.0000 
0.0000 
0.0000 
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4 Appendix 

• ecu 

1. ii~,V^,iy T ,T + 

2. ii~,V^u,d 

3. d,u,c,s 

• CC20 

4. e~,V e , 

5. e~,V e ,u,d 

• NC24 

6. z/ T ,z7 r 

7. d^v^Vp 

8. u,u,iy^,V^ 

9. /i™, r~, t + 

10. fi + ,d,d 

11. fl~ , fi + ,u,u 

12. V^V^Vr.Vr 

• NC32 

13. s, s, -u, u 

14. d, J, s, s 

15. -u, -u, c, c 

• NC19 

16. v^V^v e ,V e 

17. n~,n + ,v e ,V e 

18. It, M, Z/ e , 17 e 

19. d,d,v e ,V e 



• Mix43 

20. pT.V^v^ii' 

21. d,u,u,d 

• NC48 

22. fi~ , fi + , e~ , e + 

23. z7 M , e -,e + 

24. -u, tt, e~,e + 

25. d, J, e~,e + 

• NC64 

26. 



27. i/p, Up, Up, Up 



28. u,u,u,u 

29. d,d,d,d 

• NC25(NC24+1) 

30. b,b,iy^V^ 

31. fi~,fi + ,b,b 

• NC21(NC19+1) 

32. b,b,u e ,V e 

• NC50(NC48+2) 

33. 6, 5, e~,e + 
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